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In collisions at ultralow temperatures, molecules will pos- 
sess Feshbach resonances, foreign to ultracold atoms, whose 
virtual excited states consist of rotations of the molecules. 
We estimate the mean spacing and mean widths of these res- 
onant states, exploiting the fact the molecular collisions at 
low energy display chaotic motion. As examples, we consider 
the experimentally relevant molecules O2, OH, and PbO. The 
density of s-wave resonant states for these species is quite 
high, implying that a large number of narrow resonant states 
will exist. 



Scattering resonances are of great importance in ul- 
tracold collisions. Feshbach resonances occur when the 
energy of a pair of free atoms (or molecules) is nearly 
degenerate with that of a quasi-bound state of the pair. 
The quasi-bound state is characterized by the promotion 
of one or both atoms to an excited internal state, for 
example an excited hyperfine state Q. External mag- 
netic fields can then shift these states in or out of res- 
onance, giving the experimenter direct control over in- 
terparticle interactions In this way, Bose- Einstein 
condensates of 85 Rb have been made stable or unstable 
on command, leading to novel many-body effects ||. It 
is moreover predicted that this kind of control will be 
useful in preparing degenerate Fermi gases in "resonant 
superfluid" states with any desired interaction strength 
. Magnetic field Feshbach resonances have now been 
observed in the alkali atoms 23 Na §, 85 Rb 133 Cs 
, and 40 K |l^] . While not uncommon, these resonances 
are far from ubiquitous in ultracold alkali atoms. This 
is largely due to the relatively small number of hyperfine 
states available to form resonances. 

Ultracold molecules, by contrast, offer a far greater 
number of resonances than alkali atoms, because rota- 
tional excitations can also contribute to resonant states. 
In contrast to the two hyperfine states in an alkali atom, 
a molecule can possess a large number of energetically 
available rotational states. Considering that "typical" ro- 
tational energy splittings are of order 1-10 K, while the 
well depths of intermolccular potential energy surfaces 
(PES's) can be hundreds or thousands of K, it is clear 
that tens of rotational states may contribute, including 
their degeneracies arising from magnetic quantum num- 
bers. 

Identifying the number and properties of these reso- 
nances is important for understanding the behavior of 
a molecular gas reduced to extremely low temperatures. 



These resonances may be useful as tools for pinning down 
the details of multidimensional PES's. They may also 
imply resonant control over cold molecular collisions, in- 
cluding, possibly, control over chemical reactions |ll| . In 
addition, the scattering dynamics is chaotic, implying 
that ultracold resonances will form a new laboratory for 
studying quantum chaos |T^ | . In this Letter we estimate 
the most basic properties of these resonances, namely, 
how many we can expect, and what their widths might 
be. We are interested here in molecules produced in their 
ro-vibrational ground states by buffer-gas cooling |l3} | or 
Stark slowing jL4| techniques. However, vibrationally ex- 
cited cold molecules produced by photoassociation of cold 
atoms will also exhibit rich resonant dynamics |l5| ]. 

Classical chaos at low energies has been studied exten- 
sively in molecular scattering problems (T^| . Of greatest 
relevance to our present purposes are multiple collision 
events: it is possible that a collision deposits sufficient 
energy in internal molecular degrees of freedom (e.g., ro- 
tations and vibrations) that there is not enough transla- 
tional energy left to allow the molecules to separate. The 
molecules may therefore collide many times before finally 
shedding enough energy to separate. Multiple collisions 
are increasingly likely as the collision energy is lowered, 
and in fact are useful in quantifying the onset of chaos in 
this regime |Tsfl. 




FIG. 1. A classical trajectory for O2-O2 scattering at a 
collision energy 1 mK. See text for details. 

Our classical calculations verify that multiple-collision 
resonances persist at ultracold temperatures. We con- 
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sider scattering of 17 O2 molecules, using the singlet rigid 
rotor potential energy surface of Ref. |17| . At the collision 
energies considered, vibrational excitation is not energet- 
ically allowed. To simplify the calculation we constrain 
all the coordinates to fixed values except the intermolec- 
ular distance R and the orientations of the molecules in 
the scattering plane, as shown in the inset to Figure 1. 
The main part of the figure shows a slice through a sam- 
ple phase space trajectory for collision energy 1 mK. The 
axes are R and the angular momentum Pg\ (in units of 
h) of one of the molecules. As a guide, the heavy solid 
line indicates the allowed region of phase space. For the 
molecules to separate at this energy, the greatest angular 
momentum either molecule could have is 0.027?i. Instead, 
the molecules can exchange several h of angular momen- 
tum in each collision, and the trajectory eventually fills 
the allowed phase space. 
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FIG. 2. Evidence of chaos in quantum mechanical scat- 
tering of oxygen molecules, (a) shows the complete elastic 
scattering cross section, computed as in Ref. [19]. In (b) the 
Dyson rigidity function (points) is computed for the energy 
spacings of the 33 resonances in (a). Our results are much 
more closely related to the results of a gaussian orthogonal 
ensemble (GOE) approximation to a chaotic system than to 
a nonchaotic regular system. 

Chaotic motion in the classical realm leaves its signa- 
ture in the quantum mechanical spectrum as well. We 



have computed the quantum- mechanical elastic scatter- 
ing cross section for the \N = 0,J = l,Mj — 1) fine- 
structure state of 17 C>2 [Fig. 2(a)], following the model 
described in Ref. Jn|. As in Ref. [Q, we have neglected 
the nuclear spin of the 17 atoms, to emphasize the effect 
of from rotational motion. This model includes rotational 
levels up to N = 2, and the partial waves L = 0, 2, ...14, 
and finds 33 resonances in this energy range. Treating 
these levels as bound states of the (02)2 collision com- 
plex, we expect to find evidence of quantum chaos in 
statistical measures of the energy level distributions. A 
glance at the spectrum in Fig. 2a) suggests that there is 
no clustering of levels, but rather that they are roughly 
evenly spaced. This "rigidity" of level spacings is charac- 
teristic of the quantum mechanical spectrum of a system 
that displays classical chaos, and is usually quantified in 
terms of the spectral rigidity function A(x) of Dyson and 
Mehta |2(]] . If the number of states with energy less than 
or equal to E is denoted N(E), then A(x) represents the 
root-mean-squared deviation of N(E) from a straight line 
over an energy range x, where x is measured in units of 
the mean level spacing. Figure 2(b) shows this function, 
computed using the resonances from Fig 2(a). A spec- 
trum consisting of uncorrelated energy levels would yield 
a rigidity A(x) = x/15, as shown with the dashed line. A 
spectrum generated in the Gaussian orthogonal ensem- 
ble (GOE) approximation, representing a chaotic system, 
yields instead the solid line. The A(x) computed from 
the spectrum in Fig. 2(a) more closely resembles the 
GOE result, suggesting that this system is chaotic. A 
similar analysis has been applied to the eigenphase shifts 
in reactive scattering processes [^l] . 

Given that the classical trajectories can fill phase space 
(Fig.l), each energetically allowed state is equally likely 
to be populated during a resonant collision. This ob- 
servation has been exploited in the theory of unimolec- 
ular chemical reactions, which considers the breakup of 
a large molecule under the influence of random collisions 
with other molecules in a solution |^2j . We use the same 
idea here to estimate the total number of resonant states 
available in a given energy range. To this end we sepa- 
rate the molecule-molecule Hamiltonian into three inde- 
pendent terms: 

H = Hx + H 2 + Hint, (1) 

where Hi and H2 represent the rotational fine structure 
of each individual molecule, and H m t is an effective in- 
termolecular potential that depends only on R. Each 
Hamiltonian is represented by a set of independent en- 
ergy eigenvalues E\ , E2, and E lnt . During a resonant col- 
lision any partition of the total energy E — E\ +E2 +E[ nt 
is equally likely, provided that angular momentum is con- 
served. Estimating the number of states N(E) at or be- 
low energy E is then a simple counting exercise. The 
density of states, p(E) = dN/dE, is also easily computed. 

To illustrate that this is a reasonable estimate of the 
density of states, we compare the results to those of a 
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fully quantum mechanical calculation for 17 C>2 - 17 C>2 
scattering Table I tabulates, for the quantum me- 

chanical calculation, the number of resonances in each 
partial wave L in an energy range IK above the incident 
threshold. Each row in the table represents a separate 
calculation, characterized by the maximum number of 
molecular rotation states (-/V max ) and partial wave states 
(i max ) included in the basis set. 

Table II reports the same information, but this time 
estimated from the statistical model. Statistical fluctu- 
ations in the density of states function p(E) in the rel- 
evant energy range correspond to uncertainties of ~ 50 
% in the numbers given. Within this uncertainty the 
agreement with the quantum mechanical calculation in 
Table I is quite good, giving us further confidence that 
the statistical estimate is reasonable for this purpose. 

A clear advantage of the statistical model over the full 
calculation is that it can easily be extrapolated to include 
all allowed values of N and L consistent with conserva- 
tion of energy and angular momentum, as well as Bose 
symmetry. The result is given by the final row in Table 
II, for small values of the incident partial wave. Consider- 
ing all partial waves yields a total number of resonances 
in this IK interval of order ~ 10 3 . However, in ultra- 
cold collisions we are limited to low values of L, typically 
L = alone, barring occasional shape resonances with 
higher values of L. We therefore conclude from Table 
II that the mean density of s-wave rotational resonances 
in O2 is on the order of 10 K^ 1 . For comparison, alkali 
atoms typically exhibit on the order of 1 s-wave Feshbach 
resonance per 1 K energy interval. Thus rotational Fesh- 
bach resonances are slightly more numerous than atomic 
hyperfine Feshbach resonances, but not strikingly so, at 
least for nonpolar molecules. (In practice the numerous 
hyperfine levels in 17 02 will play a role also, but we are 
concerned here strictly with rotation.) 

This simple statistical theory can also estimate the 
widths of these resonances. Again borrowing from the 
theory of unimolecular dissociation, the resonance width 
may be approximated in the Rice-Ramsperger-Kassel- 
Marcus (RRKM) approximation as 



W{E+ 
2np 



(2) 



where p is the mean density of states we have been con- 
sidering, and W(E + ) is the number of states that per- 
mit the molecules to actually escape to infinite separa- 
tion. For the low collision energies considered here, both 
molecules must return to their ground state to separate, 
so that W(E + ) = 1. In the quantum mechanical calcu- 
lation in Fig. 2(a), there are 33 resonances in a 1 Kelvin 
interval, implying a mean width in this approximation 
of F » 5 mK. For comparison, the arithmetic mean of 
the quantum mechanical resonance widths is 8 mK. Note 
that the quantum mechanical widths span several orders 
of magnitude; their standard deviation is not given by 
the RRKM result fEl. 



The situation is dramatically different for polar 
molecules, whose long range dipole-dipole interaction po- 
tential, proportional to 1/i? 3 , holds many more inter- 
molecular bound states near threshold. To illustrate this 
point we consider two examples: OH, which is suitable 
for Stark slowing |14|; and the a(l) state of PbO, which 
is a leading candidate in which to measure the electric 
dipole moment of the electron ^5j,^6|. Both molecules 
possess, in addition to a dipole moment, a pair of closely- 
spaced states of opposite parity: a A doublet of split- 
ting A = 0.081 K in OH, and an Sl-doubling of splitting 
A = 260 /iK in PbO. Intermolecular potential curves cor- 
relating to these thresholds will hold a large number of 
dipole-bound resonant states. Details of the fine struc- 
ture of OH is given in J||], and that of PbO in |f]. 

For the OH-OH interaction an approximate PES ex- 
ists Jm!. We include the long-range part of this inter- 
action, which has a minimum at R w 6 a.u. that arises 
from hydrogen bonding forces. For the sake of our sta- 
tistical model, we approximate this surface with a sin- 
gle curve by fixing the two OH molecules in their lin- 
ear geometry. The intermolecular bound state energies 
-Bint are then simply the bound states of this potential, 
for each of the partial waves considered. We ignore for 
our present purposes an even deeper part of the PES 
corresponding to chemical binding of hydrogen peroxide, 
H202- For PbO-PbO the PES is entirely unknown; how- 
ever, since the bound states of interest depend primarily 
on the long-range potential, we estimate the PbO-PbO 
PES by multiplying that of OH-OH by the appropriate 
ratio of squared dipole moments, [dpho/^on] 2 ~ 7.7. Ta- 
ble III tabulates the number of s-wave resonances to be 
expected for each species. These numbers refer to the 
important energy range 2A above threshold, where the 
number of resonances is greatest. Expressed as a mean 
density-of-states over this range, the densities are 240 
K _1 for OH and 2 x 10 5 K" 1 for PbO, significantly larger 
than for the nonpolar O2. 

In general, the resonances shift when an external elec- 
tric field is applied, so that one after the other can be 
brought into resonance with the essentially zero-energy 
molecules. Given a mean density of states w 10 2 K _1 
for OH, and given the OH dipole moment (=1.67 D), the 
mean spacing of these resonances as a function of field 
would be of order ~ 250 V/cm. This means that, in 
contrast to the magnetic field resonances in alkali atoms, 
there will be many available resonances in ultracold polar 
molecule collisions. 

We conclude by sketching the implications of this work 
for the theory of ultracold collisions. As is well-known, 
ultracold collision theory requires a certain level of exper- 
imental input to fine-tune the scattering Hamiltonian. 
It is not yet clear how much or what kind of informa- 
tion is required to perform this fine-tuning for molecules 
p9| . Knowing roughly how many resonances to expect 
can greatly simplify theory, by providing an estimate of 
how many parameters are strictly necessary to reproduce 
the data. With this information theory could then turn 
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around and predict how these resonances evolve with ap- 
plied external electric and magnetic fields. 

A second, equally important implication concerns ex- 
periments at temperatures where kT is large compared 
to the mean spacing between resonances. In this case in- 
dividual resonances are not resolved, and only the overall 
modulation of the spectrum matters. It is then likely that 
semiclassical methods should prove insightful. 
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TABLE I. Number of resonances in a IK energy interval 
above threshold for cold collisions of 17 02 molecules in their 
J = M.j = +1 state, neglecting nuclear spin. These results 
are determined from a quantum mechanical calculation. Var- 
ious levels of completeness of this calculation are specified by 
the maximum value of rotational quantum number A max and 
maximum partial wave L max . For each such calculation, the 
table gives the number of resonances found separately in each 
incident partial wave L. 

AT max L max L = L = 2 L = 4 L = 6 L = 8 L = 10 
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TABLE II. Number of resonances for 1 O2 cold collisions, 
as in Table I. Here, however, the resonances are estimated us- 
ing the statistical counting model described in the text. In ad- 
dition, the final row extrapolates to the limit where A max = 8 
and I/max = 36, thus including all possible energetically al- 
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TABLE III. Summary of s-wave resonances for two po- 
lar molecules. pA is the average density of states consistent 
with s-wave initial conditions, over an energy range 2A above 
threshold, where A refers to the A (n)-doublet of OH (PbO). 
Using the RRKM result, this density translates into the mean 
energy width Pa- 
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A(K) 


PA (K- 1 ) 


f A (A»K) 
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